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We consider the problem of estimating convex boundaries from 
blurred and noisy observations. In our model, the convolution of an 
intensity function / is observed with additive Gaussian white noise. 
The function / is assumed to have convex support G whose bound- 
ary is to be recovered. Rather than directly estimating the intensity 
function, we develop a procedure which is based on estimating the 
support function of the set G. This approach is closely related to the 
method of geometric hyperplane probing, a well-known technique in 
computer vision applications. We establish bounds that reveal how 
the estimation accuracy depends on the ill-posedness of the convo- 
lution operator and the behavior of the intensity function near the 
boundary. 

1. Introduction. We consider the problem of recovering convex bound- 
aries of an image from blurred and noisy observations. Following Hall and 
Koch [18], we concentrate on the so-called continuous image model (which 
is explained in what follows). The original image is given by an intensity 
function f, supported on a closed convex set G in the d-dimensional Eu- 
clidean space M. d . In many instances, only a degraded version of the original 
image is available, where the typical reasons for degradation are blurring 
and noise; we refer to [4], pages 51-53, for a detailed discussion of different 
sources of degradation in imaging. Blurring is modeled using a convolution 
operation, K * /, in which the original image / is subjected to the effects of 
a point spread function K. The effects of noise degradation can then be ad- 
equately captured by adding a stochastic component to the blurred image. 
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In this setup, the standard problem is to recover the original image / from 
its degraded version. 

The outlined inverse problem of image deconvolution is the subject of a 
considerable literature in image processing and statistics. In this paper, how- 
ever, we focus on directly recovering the boundary of the image, rather than 
estimating the image itself /. We note that in virtually all image processing 
applications, edge detection is one of the standard preliminary analysis steps; 
this fact serves to motivate our paper. Although many practically important 
edge detection algorithms for degraded images have been proposed in the 
image processing literature, their theoretical properties are rarely analyzed. 
Our goal in this paper is to study the theoretical properties of a particular 
algorithm which uses geometric probing to estimate the boundary of the 
image. 

Our paper is closely related to two strands of research. The first focuses on 
nonparametric signal and image deconvolution, which has been extensively 
studied in statistics; see, for example, [10, 11, 14, 17, 18, 22] and references 
therein. In this literature, the function / is typically assumed to be smooth, 
even though its smoothness may be unknown. The second stream of research 
considers the problem of recovering the boundary of a multidimensional 
image from direct observations; see, for example, [8, 19, 20, 21, 22, 25, 26, 31] 
and references therein. Recovering boundaries in models that involve indirect 
observations has been discussed recently in [7] and [15]. 

The main contributions of this paper are the following. We propose a 
method for estimating the boundary of the support G of an image /, where 
information on the image is only available in the form of blurred and noisy 
observations. Our approach focuses on direct recovery of the support func- 
tion of the set G; the proposed estimation scheme is akin to the method 
of geometric hyperplane probing which is very common in computer vision 
applications (see, e.g., [30] and [24]). A similar estimation method has been 
proposed in [15] for reconstructing convex shapes from noisy observations of 
their moments. In terms of accuracy of the proposed estimation procedure, 
we first study the behavior of an estimator of the linear probe functional 
(see Theorem 1), and subsequently establish upper bounds on accuracy of 
the support function estimation algorithm (see Theorem 2). We next extend 
this pointwise result by deriving upper bounds on the global estimation ac- 
curacy, as measured by the Hausdorff distance between the boundary and 
its support-function-based estimate (see Theorem 3). 

It is worth noting that an alternative approach to the problem studied 
in this paper is to first deconvolve the function /, and subsequently infer 
the support boundary. However, in our setup the assumptions imposed on 
the intensity function / are quite weak. In particular, we assume that / is a 
square integrable bounded function with convex support, with some restric- 
tions imposed on its behavior near the boundary. Under such assumptions, 
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the accuracy of any deconvolution estimator can be arbitrarily bad. Hence, 
in the absence of additional information on the regularity of /, "direct" 
estimation of the boundary can be advantageous. 

The paper is organized in the following way. In Section 2 we formulate 
the estimation problem, and introduce some preliminaries and notation. Sec- 
tion 3 details our estimation procedure, whose accuracy is analyzed in Sec- 
tion 4. Section 5 contains some discussion, and proofs are all relegated to 
Section 6. 

2. Problem formulation. 

The model and observation process. The image is modeled as an un- 
known positive function / G L2(R d ), the so-called intensity function, sup- 
ported on the convex compact set G with nonempty interior in R d . Suppose 
that we can observe the process Y(x), given by 

(1) dY(x) = {Kf)(x)dx + edW{x), xeR d , 

where K:T>{K) \— > 1Z(K) is a linear transformation with domain T>(K) C 
L 2 (R d ), and range K(K) C L 2 (R d ), < e < 1, and (W(x) : x £ R d ) is a stan- 
dard d-dimensional Brownian motion. The observation scheme (1) implies 
that, for any function q £ L 2 (R d ), we can observe (q,Kf) with added zero 
mean Gaussian noise having variance e 2 ||g|| 2 . Here, (•,•) and || • || denote 
the inner product and the corresponding norm in L 2 (R d ) with respect to 
Lebesgue measure. In what follows, we restrict attention to the case where 
K is the convolution operator 

(Kf)(x)= I K(y)f(x-y)dy, 

(2) Jr " 

x£R d j£L 2 (R d ),K&L 1 (R d ). 

The function K : R d \— > R d which models the blurring process is referred to 
as the point spread function in image analysis (cf., e.g., [4], pages 51-53). 
Following Hall and Koch [18], we refer to (l)-(2) as the continuous image 
model. Our goal is to estimate the boundary of G, using observations from 
(1)"(2). 

Problem geometry and notation. We introduce some notation that will 
be used in the sequel. The Euclidean norm of the d- vector x = (x\, . . . ,Xd)' 
is denoted by \x\. Let B d = {x:|x| < 1} denote the centered unit ball in 
R d ; its surface is the unit d-sphere 3 d " 1 = {x : \x\ = 1}. For a unit vector 
u G 5 d_1 , we denote by u 1 - the (d— 1) -dimensional subspace orthogonal to u. 
The centered cube with faces parallel to the coordinate hyperplanes will be 
designated by E = {x : \xi\ < 1, 1 < i < d}. In what follows we will consider 
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some simple affine transformations in M. d . In particular, for u£ S d 1 , let A u 
denote the d x d rotation matrix which maps (1,0,..., 0)' to the unit vector 
u G S d ~ l . For fixed r = (u,r) G T := S^" 1 x [0, 1], we define 

(3) E T := A U E + (1 + r)u. 

Thus, E T is obtained by rotating E according to A u , followed by a translation 
x i— ► x + (1 + r)u. Consequently, -EV is the cube centered at (1 + r)u with 
sides of length 2 and faces parallel to the coordinate hyperplanes of (u,u ). 

We assume that G, the support of /, is a cf-dimensional convex body, 
G C The following concept of a support function plays a key role in our 
estimation procedure. 

Definition 1. For u G S d_1 , the support function h(u) of the set G is 
defined by 

h(u) := sup{a/ti : x G £?}. 

The supporting hyperplane to G with outward normal u G S^" 1 is given 
by {x : x'u = h(u)}; the support function h(u) at the unit vector u gives the 
signed distance from o to the supporting hyperplane. Throughout the paper 
we assume that the origin, o, belongs to the interior of G. This assumption is 
not restrictive; if it holds, the support function h(u) gives the actual distance 
from o to the supporting hyperplane. 

The key observation that will be exploited in what follows is that there is 
a one-to-one correspondence between a closed convex set G and its support 
function, namely, 

(4) G = {x:x'u<h(u),VueS d - 1 }. 
For fixed u G S^" 1 and small n > 0, we denote 

(5) H u (n) := {x : h{u) -n<x'u< h(u)}. 

In words, H u {rf) is the set formed by the intersection of G with the half-space 
{x : x'u > h(u) — 7]}. The volume of the set H u {rf) characterizes "massiveness" 
of G in the direction u in a small vicinity of the support value. It is interesting 
to note that behavior of the volume of H u (n) as n — > is related to the tail 
behavior of the Fourier transform of Ig(') 111 the direction u G S^ -1 (see, 
e.g., the survey paper by Brandolini, Rigoli and Travaglini [5] and references 
therein). We remark also that H u {rj) = H G because G C B d ; this 

will be used repeatedly in what follows. 
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3. The proposed estimation scheme. Our approach to recovering the 
boundary of the set G is based on pointwise estimation of its support func- 
tion h{-). Before developing the details of the proposed approach, let us first 
highlight the main idea that underlies it. Fix a direction u, and imagine that 
a (i-dimensional cube is placed outside the unit ball in R d so that one of its 
faces is tangent to the ball in direction u. This will be referred to as the cube 
probe. Recall that G C B , so by "sliding" this cube toward the origin along 
the axis u, the cube eventually intersects with the boundary of G. In the 
absence of blurring and noise (i.e., K = I and e = 0), the cube probe will 
"hit" the boundary when it is exactly at distance h{u) from the origin. As 
the probe penetrates G, it accumulates "mass," where the rate of this mass 
accumulation is determined by the following: (i) the intensity function (/); 
(ii) the blurring incurred by the convolution operator (K); and (iii) the noise 
level (e). The main idea is to estimate the distance where the accumulated 
mass begins to differ significantly from zero, indicating the location of a 
boundary. 

The probe functional. For fixed r = (it, r) G T, where T := S"^ -1 x [0, 1], 
let E T be given by (3). We define the probe functional by 

£ f (r):= f f(x)dx = (f,l Er ), T = (u,r)eT, 

where 1b stands for the indicator function of a set B G M. d . The value of 
the linear functional £j{t) is nothing but the mass accumulated by the cube 
probe when it is at the location r = (u, r) G T. The important property of 
the probe functional which underlies our construction is the following. For 
any fixed u G S^ -1 , if r G [h(u),l], then £t(u,r) =0; as r decreases from 
h(u) to 0, £f(u,r) grows monotonically, because / is positive. Although we 
suppose that / is everywhere positive, it is seen from the proofs below that 
it suffices to assume that / is positive near its support boundary, that is, 
for u G S d-1 , there exists A > such that f(x) > for all x G H U (A). 

Estimation of the probe functional. The first step of our construction is 
to estimate the probe functional £f(r) for given r = (u,r) G T. Fix 5 > 0, 
and let p$ denote a real- valued infinitely differentiable function p$ G C°° (M 1 ) , 
equal to 1 on the interval [— 1 + 5, 1 — 5] and outside [—1, 1]. Define 

d 

(6) <P8(x):=Y[ps(xj), x = (x 1 ,...,x d )'. 

3=1 

Clearly, the support of ipg is the cube E = {x : \xi\ < 1, 1 < i < n}. For fixed 
r = (u,r) G T, we define the shifted and rotated version of ip$ to be 



(7) 



Vt,s{x) := ips(A u x - (1 + r)u). 
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By construction, the support of (p Tj g is the shifted and rotated cube, E T , 
defined in (3). 

Suppose that cp T; $ is in the range of K* , TZ(K*), where K* is the adjoint 
operator to K. Then, according to the linear functional strategy (cf. [1] and 
[13]), if V(K) is dense in L 2 (R d ), then there exists Vv,«5 G L 2 (M d ) such that 
ip Tt S = K*ip Tt s, and 

(8) (f,<PT,s) = (Kf,i/>T,8) = (f,K*il>T,8) V/GD(K). 

Now, for any function q £ Li(M d )nL2(K d ), let q denote its Fourier transform 

q(u):= q(x)e w ' x dx, u = (ui, . . . ,u d )'. 

It is easily verified that, for the convolution operator K, the range of the 
adjoint is given by 

■Jl(K*) = ^q: J d \q{uj)\ 2 \K{uj)\~ 2 aluj < ooj, 
and Vv,<5 m (8) is given by 

(9) M^i/,^^- 

The following result states that, under natural assumptions on K, ip Tt s 
belongs to 1Z{K*) for any r G T and all 5 > 0. 

Assumption 1. There exist constants L > and /3 > such that 

|K(w)| >L(i + |w| 2 )-' 3 / 2 Vwet d . 

Lemma 1. Lei Assumption 1 /ioW. Then, (p Tt $ E 1Z(K*) for any t G T 
and aZZ <5 > 0, and in addition, 

(10) ||^|| a <GiL- a rw 

where C\ = C\(d,[3) depends on d and (5 only. 

We note that Assumption 1 guarantees the identifiability of / from the 
observations (l)-(2), and states that the tails of the Fourier transform of the 
point spread function K cannot decrease to zero faster than the indicated 
polynomial rate. This assumption is quite standard in deconvolution prob- 
lems, and corresponds to what is known as a moderately ill-posed problem. 
We note that the severely ill-posed case, where the tails of K are allowed to 
decrease at an exponential rate, can be also treated using our methodology. 
This setup is, however, beyond the scope of our paper. 
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We are now ready to define the estimate of the probe functional £f(r) 
based on the observations (1). For fixed r = (u,r) G T and 5 > 0, we define 

(11) £{t;6):= f ^ TyS (x)dY(x), 

where ip T 5 is given in (9). We note that, in view of Lemma 1, the estimator 
£(t;5) is well defined for all r G T and 5 > 0. In what follows, we assume 
that 

(12) sup I /(x) I <M<oo. 

x&G 

The next statement establishes an upper bound on the accuracy of the probe 
functional estimator £(t;5) for a fixed direction u G S"^" 1 uniformly over 
re [0,1]. 

Theorem 1. Let Assumption 1 and (12) hold, and suppose that e is 
sufficiently small. Define 

1/(0+1/2) 

(13) 




and £*(t) :=£(t;5*). Then for any fixed u G S d 1 , we have 

r 1V2 / r^r\ V(/?+i/2) 

(14) |E sup |4(r)-^(T)| 2 | ^^K^fM^-^ln-j 

It is important to emphasize that the upper bound of Theorem 1 holds for 
all square integrable functions / supported on G C B d and satisfying (12). If 
further assumptions on / are introduced, for instance, if / is assumed to be 
smooth, then the rate of convergence in estimating £f(r) may be improved; 
see further discussion in Section 5. 

The support function estimator. Based on the estimate £*(t) = &*(u,r) 
of the probe functional £f(r), we define the estimator of the support function 
h(-) at a fixed direction u G S^" 1 in the following way: 

(15) h = h(u) :=max{r G [0, 1] :£*(u,r) > 9}, 

where £*(u,r) = £*{t) is given in Theorem 1, and 6 > is a threshold to be 
chosen in the sequel. If £*(u, r) < 6 for all r G [0, 1], we set h(u) = 0. 

We note that, for 9 sufficiently small, the estimator h(u) is well defined. 
Indeed, for any fixed u G S^ 1 , £*(r) =£*(u,r), considered as a function of 
r, is a separable Gaussian process on [0, 1] that has continuous sample paths 
with probability 1. The last fact follows from the bounds in the proof of 



8 



A. GOLDENSHLUGER AND A. ZEEVI 



Lemma 2 in Section 6, continuity properties of ip T g as a function of r and 
the well-known criteria for continuity of sample paths of Gaussian processes 
(see, e.g., [12], pages 193-194). In particular, this implies that £*(u,h) = 9 
with probability 1. 

4. Accuracy of the support function estimator. 

The class of intensity functions. To analyze the accuracy of the estimator 
h(u), we introduce the following class of functions /. 

Definition 2. We say that a positive function / G L2(M. d ) with convex 
support G C B d , o G int(G), belongs to the class T u (ot), with a > 1, if (12) 
holds and for a fixed direction u G 5 rf_1 , there exist Q > 0, A > such that 

(16) / f(x)dx>Qrj a Vr?G(0,A), 
where H u (r]) is defined in (5). 

Several comments on the above definition are in order. The integral on 
the left-hand side of (16) quantifies the "massiveness" of the image / in the 
direction u. In view of our definition of the probe functional, this integral is 
nothing but £f(u,h(u) — r/). Thus, (16) can be equivalently rewritten as 

(17) £ f (r) =£ f (u,r) >Q\h(u) - r\ a Vr G (h(u) — A, h(u)). 

It is important to realize that the class J- U {a) is defined for a fixed direc- 
tion u G S and therefore, the constants Q, A and a may depend on u. 
Relation (17) imposes restrictions both on the geometry of G in the vicinity 
of the support value h(u) in direction u G S d ~ l , and on the behavior of / 
near the boundary of G. The probe functional £f(u,r) quantifies the rate at 
which the cube probe E T = Er U} r) accumulates mass as it penetrates the set 
G in a fixed direction u (as r decreases from 1 to 0). The rate of increase 
in accumulated mass depends on the behavior of the intensity function / 
near the boundary, and the local curvature of the boundary itself. For ex- 
ample, if / has a discontinuity jump on the boundary, that is, / = lcf for 
/ > c > 0, the probe functional £f(u,r) behaves roughly like the volume of 
the set H u (h(u) — r) = E T n G. For simple geometrical objects, the order 
of this volume can be easily derived. In general, it turns out that the or- 
der of this volume is essentially determined by the (Gauss) curvature of the 
boundary (see, e.g., [5, 6]). We note that the parameter a indexing the class 
J-" u (a) satisfies a > 1 whenever G is a convex set. The next examples illus- 
trate how the curvature of the boundary of G and the behavior of / near 
the boundary determine the index a. In what follows we use the term sharp 
boundary when the intensity function / is discontinuous on the boundary of 
the support set. 
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Example 1. Assume that d = 2 and consider the case of a sharp bound- 
ary. Let G be a convex polygon, and for simplicity, suppose that f = 1q- 
Then for any direction which is not perpendicular to the side of the polygon, 
we have a = 2. In this case the corresponding supporting line intersects the 
boundary of G in a vertex of the polygon, and the constant Q depends in 
an obvious way on the angle between two adjacent sides of the vertex. If 
the direction uGS 1 is perpendicular to the side of the polygon, then the 
supporting line contains that side; in this case a = 1. 

Example 2. We again consider the case of d = 2, and first assume that 
the boundary is sharp. If G is a circle or an ellipse, then a = 3/2 for any 
direction. It turns out that this case is rather general. Let x u denote the point 
on the boundary dG where the support value in direction u is attained. If dG 
has nonzero curvature at x u , then a = 3/2 for that direction. Now assume 
that the boundary is nonsharp. Specifically, without loss of generality, we 
let u= (0,1) E S 1 , and suppose that dG can be represented as X2 = g(x\) := 
— x\ + c, with c > 0, in the vicinity of the origin. This corresponds to the case 
of nonzero boundary curvature in direction u. Let f{x\,x<i) > \x2 — g(xi)\" / 
for some 7 > (7 = corresponds to the sharp boundary case). Then it is 
easily verified that = 7 + 3/2. If the curvature of dG vanishes at x u , the 
exponent a differs from 7 + 3/2. For instance, if dG is represented as the 
graph of the function X2 = —x 2 ^ + c, q > 1, then a = 7 + 1 + (2q)~ 1 . We 
note, however, that if the boundary is smooth, then the set of points on the 
boundary where the curvature vanishes has zero Lebesgue measure. 

Example 3. In M. d we consider the case of a sharp boundary. If G is an 
ellipsoid in ]R rf , then a = (d + l)/2 for any direction u G Similarly to 

the previous example, if for some u £ S d ~ l the Gauss curvature of dG does 
not vanish at x u , then a = (d+l)/2 as well. If G is a parallelotope, a takes 
values in the set {1,2, ... ,d} depending on the direction u. In particular, 
if u is perpendicular to a face of the parallelotope, then a = 1. Without 
loss of generality, consider now the situation when u = (0, ... ,0, 1) € S^" 1 
and dG can be represented as x^ = g(x\, . . . , x^-i) ■= —x^ gi ■ ■ ■ xf^ 1 + c for 
some qi > 1, i = 1, . . . , d— 1, and c> 0. Then it is not difficult to verify that 

« = i + E J ti 1 (2g J )- 1 . 

Accuracy bounds. We are now ready to establish upper bounds on the 
accuracy of the proposed estimation procedure. 

Theorem 2. Let Assumption 1 hold. Fix an arbitrary u G S , and 
let h*{u) be the estimator of the support function given by (15), with 9 = 9* 
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taken to be 

(18) 9*:= ( jM^ 2 yJc 3 \n- J , C 3 > 0. 
Then, for sufficiently small e, 

( I =- i l/(a(/9+l/2)) 

sup {E|^(u)-h(«)| 2 }V 2 <(74Q-V«JM^-V 2 £Jl n - I 

where C3 and C4 may depend on d, (3 and A on/y. 

We note that since the blurring kernel K is assumed to be known, the 
dependence of the above threshold 9* on L and (3, which characterize the 
behavior of K, is not restrictive. 

Theorem 2 gives an upper bound on the pointwise risk of the proposed 
estimator. Based on the estimate h*, we can define the global estimator of 
the set G as 

(19) G*:={x:x'u<h*(u),VueS d - 1 }. 

We note that is a convex set, by construction, because it is given by the 
intersection of half-spaces formed by the estimated supporting hyperplanes. 
In order to measure global accuracy of G* as an estimator of G, it is natural 
to use the following family of L p -metrics on the class of all convex bodies 
in M. d . If G\ and G2 are convex bodies with support functions h Gl and hc 2 i 
respectively, then the L p -metric is defined as 

A P (G 1 ,G 2 ) := y \h Gl (u) - h G2 (u)\ p dXa-xiu)] 11 \ 

where A^_i is the spherical Lebesgue measure on S^ 1 . It is remarkable that 
the Aoo-metric is nothing but the usual Hausdorff distance between the sets 
(see, e.g., [29], Section 1.8). We note also that the metric A p may be viewed 
as the general L p -metric of [3] specialized to convex sets; the cited paper 
argues that such a metric has attractive properties for image processing 
applications. 

To state the upper bound on the risk of G* we first introduce the following 
"global" version of the functional class J- U (a). 

Definition 3. We say that a positive function / with convex support 
G<ZB d , oG int(G), belongs to the class J~(a), if (12) holds and for all 
u E there exist Q > 0, A > and a > 1 such that 



1 < P < 00, 



where H u (r]) is defined in (5). 



/ f(x)dx>Q V a V 7/6(0, A), 

JHJri) 
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In contrast to Definition 2, here we require that (16) holds for all directions 
u G S^ 1 with the same constants Q, a and A. This class may be adequate, 
for example, for purposes of describing indicator functions of convex sets 
with smooth boundaries having everywhere nonvanishing Gauss curvature. 

It follows immediately from Theorem 2 that 



sup {EA^(a,G)} 1/2 <C 5 Q" 




l/(a(/M-l/2)) 



Our next result concerns global estimation accuracy as measured by the 
HausdorfF distance. 

Theorem 3. Let G7* be given by (19). Then, under the conditions of 
Theorem 2, 

( I =-<) l/W+l/2)) 

sup {EAL(G„G)} 1 / 2 < C % Q- l ' a \ M^L-Wlni 
fGF{a) { V 6 J 

5. Discussion. We now turn to a few brief comments on the main results 
of this paper. 

1. The proposed algorithm is based on direct estimation of the edge by 
means of the geometric hyperplane probing scheme. As our results show, the 
accuracy of this method is determined by the ill-posedness of the convolution 
operator K (given by the value of /3) and the "massiveness" of the image 
in the u-direction (near the boundary) as quantified by the index a. The 
latter is determined by the local behavior of the intensity function / near 
the boundary and the local geometrical properties of the support G in the 
vicinity of the estimated support value; these two factors determine the 
index a in the manner illustrated in Examples 1-3 in Section 4. 

2. The functional class J- U (a) involves only weak assumptions on the be- 
havior of the intensity function. Specifically, / £ J- U {ct) is only assumed to be 
bounded and square integrable with convex support, satisfying certain "mas- 
siveness" conditions near the boundary. We would like to emphasize that, 
under these conditions uniform bounds on the accuracy of a deconvolution- 
based estimator can be arbitrarily bad. This is due to the fact that the 
functional class J- u {o) is too broad. Hence, our approach, which is built on 
"directly" estimating the boundary, can be applied in instances where the 
deconvolution-based estimator is not appropriate. 

3. The basic ingredient in our construction is the estimation of the probe 
functional corresponding to the hyperplane probing scheme. Using Theorem 
1 and results in [9] , we can show that our estimator of the probe functional is 
rate optimal (up to a logarithmic term). The proposed boundary estimator 
inherits this convergence rate. 
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4. If additional regularity conditions on / are imposed, then one can im- 
prove the rates given in the main results of this paper. In particular, there 
is a close connection between the setup and problem formulation in our pa- 
per, and that of estimating a change-point from indirect observations in the 
one-dimensional case (see, e.g., [27]). The recent paper by Goldenshluger, 
Tsybakov and Zeevi [16] studies the problem of change-point estimation 
in a function / in the white noise convolution setup. In that paper it is 
shown that minimax rates of convergence are determined by smoothness 
of / away from the change-point, and by the ill-posedness of the convolu- 
tion operator. Specifically, further smoothness of the function / away from 
the change-point results in better accuracy in estimating the change-point. 
(This stands in contrast to the case of direct observations.) This suggests 
that the optimal rate of convergence in estimating convex boundaries from 
indirect observations depends on the smoothness of the intensity function 
/. Establishing this result rigorously remains an open problem. 

6. Proofs. 



PROOF of Lemma 1. We assume that 5 > is fixed, and note that, by 
(7), 



(20) 



ip Tt s(x)e tu x dx 
ips(A u x - (1 + r)u)e* ( 



' dx 



dy 



= exp{z(l + r)u'A u Lo} 
= exp{i(l + r)v! A u u}tp s (A u u). 

It follows from (6) that (p$(oS) = Ylj =1 ps(ujj) for all u = (u±, . . . ,uJd)' G 
The function p$(-) is rapidly decreasing, that is, |ps(A)| < Cfc(l + |A|)~ fc for 
all k. Using this fact, and taking into account Assumption 1 and (20), we 
conclude that Cp TjS (-)/K{—) G n L 2 (R d ) so that ip Tl s(-) is well defined 

in (9). This also implies that (p Tt g G 1Z(K*). 

Now we prove (10). Straightforward algebra shows that 

(ps(A u u) I 2 



(21) 



<L" 



L 



-2 



dw 

K{-u>) ■ 

\(p & {A u uj)\ 2 {l + \u>\ 2 YdLO 

\(p s {uj)\ 2 {l + \uj\ 2 fdjjj 



niwK)i 2 (i+^i 2 )^ 

i=l 

|^(A)| 2 |A| 2 ^A, 
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where c\ = c\{j3,d) is the constant depending on d and only. Thus, it 
is sufficient to bound the last integral in (21). To this end, we recall the 
standard way to construct the function ps G C°°(]R 1 ) with aforementioned 
properties (see, e.g., [23] and [28]). Let j(x) = e~ 1 ' I^ 1-2 ^ . Then, on the 
interval [—1,-1 + 5], where ps climbs from zero to one, set 

ps(x) = j F_i(^r) d v> *e [-1,-1 + 5], 
where C2 is an absolute constant. First assume that (3 is an integer; then 



and therefore 

r ra\ n r£ r in in/.t-i-IX 2 

dx 



[\pf\ X )\ 2 d X = ^[ 

JR1 Z P 7R1 



6 



C 2_f ^-^(xtfdx 



< 



5 2 /3- 
C3 



§2/3-1 ' 

where C3 = cs(P) depends on (3 only. Combining this with (21), we obtain 
(10) for integer j3. For general (5, the result follows from the standard inter- 
polation inequalities for Sobolev spaces (cf. [2], page 127). □ 

We now provide an auxiliary lemma that will be used repeatedly in the 
proofs below. Define 

(22) X( T )=X(u,r):=e [ V+<5, (x) dW(x), re [0,1], 

JR d 

where r = (u,r) £ T, ijj Tt s t is given by (9) and 5* is defined in (13). In the 
sequel we will treat X{r) as a random process over the index set T = 5 ,rf_1 x 
[0,1]. When u £ S d ~ l is fixed, we have the random process {X(u,r) :r G 
[0, 1]}. This distinction will always be clear from the context. 
Obviously, {X(t)} is a zero mean Gaussian process. Let 

a|:=supE|X(r)| 2 , 

rGT 

and note that, in view of Lemma 1 and (13), we have that 

2 2|| ; 1 1 2 

a x = supe ||V+<ul 

reT 
CiE 2 

(23) - 2/3 _! 



£ 2 N l/(/3+l/2) 
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Lemma 2. Let Assumption 1 and condition (12) hold. 
(i) For any fixed u £ S^" 1 and for all 6 > 2a x, we have that 

p( sup \X(u,r)\ >e\< ClI f e exp(- 
Ue[o,i] J a x < 



E[0 

where C\ = C\{d,(3) depends on d and (5 only. 
(ii) For all > 2o"x , we have that 



2a| 



A' 



exp 



24 



p|sup|X(r)| > (9 J < C : 

where C2 may depend on d, (3, M and L only. 

Proof, (i) Define the semi- metric 

p(p,r) :=(E\X(u,p)-X(u,r)\ 2 ) 1/2 . 

To prove the inequality in the lemma, we verify the conditions of Proposi- 
tion A. 2. 7 in [32]. Specifically, we need to evaluate the minimal number of 
balls of radius v with respect to the semi-metric /i covering the index set 
[0,1]. We have that 



p 2 (p,r) 



= £ 

< £ x 

(a) , 

< e 1 

(b) 2, 

< £ 2 \p-r\ 2 



^ 2| |2 

< e \p — r\ 



(c) 

< ci I/O — r\ 



«,p),<5* V^i^r),^ H2 



A(-u,) 

^0,p),5, M - £(u,r),A, ( w )| 2 (! + \u\ 2 f duJ 



\ip Sf (A u u;)\ 2 (l + MY |1 - j(p-r)WA v uf ^ 



\<ps* M| 2 (i + M 2 yVw| 2 dw 



= ciM 2 \p-r\ 2 , 

where (a) follows from (20); (b) follows from a change of variable noting 
that A u is a rotation matrix, and the standard bound [1 — exp{ix}| < 
(c) follows from the same argument as in the proof of Lemma 1; and the 
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last equality follows by definition of 5* given in (13). Therefore, the min- 
imal number of balls of radius u, with respect to the semi-metric p, that 
covers the interval [0,1] does not exceed N(v; p, [0, 1]) = C2Mu~ 1 . Applying 
Proposition A. 2. 7 from [32] with eo = &x, we obtain that, for some constant 
c 3 >0, 

f 1 c 3 M9 f 9 2 

F{ sup \X(u,r)\ > 9 } < exp<^ -— =- 

Ue[o,l] > a x >- la x 

provided that 9 > 2a x ■ This concludes the proof of (i) . 
(ii) Now we consider 

p(r,t) = (E\X(r)-X(t)\ 2 ) 1/2 , 

where t = (u, p), t = (v,r) and r,(6T = 5' d_1 x [0,1]. Then, similarly to the 
above considerations, we have 



p 2 (r,t) < e 2 \$ rA {uj)-$ tA {Lj)\ 2 {l + \uj\ 2 f<ko 

< Ci e 2 l I [^(A^l'le^+^'^-e^+^'^l^l + lu;! 2 )^^ 
+ / [^(^^-^(^^l^l + lwl 2 )^^ 



=: c 4 e 2 {J 1 + J 2 }- 
Straightforward algebra shows that 

■h<c 5 [ \(p St {A u uj)\ 2 \{l + r)u' A u uj - {I + p)v' A v lo\ 2 {I + \uo\ 2 f duj 

JR d 

= c 5 (l + p) 2 

x / \$ 5f (A u uj)\ 2 \jA' u {u-v)+J(A u -A v )'v- r -f- P -JA' v v 

JR d I 1 + p 

(24) x (1 + \uj\ 2 ) p dw 

<c 6 / \(p s ^A u u)\ 2 \u\ 2 {\u-v\ 2 + {r- P ) 2 }{l + \u\ 2 fdu) 



<c 7 \t-r\ 2 \(ps t (A u u)\ 2 (l + \u 



dw 



<c s \t-r\ 2 (L5f +1 r\ 



where (a) follows from the fact that A u and A v are orthogonal matrices, 
hence, from a singular value decomposition, one has that the spectral norm 
of the difference between A u and A v is 0(\u — v\) when u,v € S^" 1 are close; 
and (b) is obtained as in (i). 

In order to bound J2 from above, we first assume that /3 is an inte- 
ger. We introduce the standard multi-index notation: D k = D^ 1 ■ ■ ■ D^ d = 
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«9l fc l /dx kl ■ ■ ■ dx k d d is the differential operator of order | k\ = k\ + • • • + kd, 
k = {k\ , . . . , kd) , and for x = (x\ , . . . , Xd) , we write x k = x kl ■ ■ ■ x k d d . Further, 
we note that <fs(') 1S infinitely differ entiable; hence, there exists < ( < 1 
such that 

(25) (ps, {A u lo) = (p s „ (A v oj) + {A v u - A u lo)' V<ps» (A v co + (A u u;), 
where V := (D\, . . . ,D P ) is the gradient vector. Using (25), we can write 

(26) J 2 <c 8 \u-v\ 2 f \V<psA(A v + (A u )u})\ 2 (l + \u>\ 2 f +1 du;, 

JR d 

and our current goal is to bound the last integral. For u and v close to each 
other, the matrix I + ( > A' V A U is nonsingular, thus, using a change-of- variables 
in the above integral, we see that it suffices to bound from above the integral 
J3 := / K d I Vips, ( w )| 2 (l + |w| 2 )^ +1 du>. To this end, we note that 

J 3 <c 9 £ / \D k {(-ix)<psAx)}\ 2 dx 
\k\^+i Jud 



< cio J \Pg dx < cn<5, 

where the last two inequalities follow from the definition of ips„ (and ps„ ) and 
from the same argument as in the proof of Lemma 1. Combining this with 
(24) and (26), and substituting (13) for 5*, we obtain /i 2 (t,r) < cu\t — t\ 2 . 
Thus, the minimal number of balls of the radius u in the semi-metric \i cov- 
ering the index set T = S"^ -1 x [0,1] does not exceed N(u;^i,T) < cizv~ d . 
Therefore, applying Proposition A. 2. 7 from [32], we come to the result (ii). 
The statement of the lemma for noninteger [3 follows from standard inter- 
polation inequalities for Sobolev spaces (cf. [2], page 127). □ 

Proof of Theorem 1. It follows from (1), (11) and (8) that 
£(t; 5) = (Vv,<5, Kf) + e f ^s(x) dW{x) 
= (f, <Pt,s) +£ i>r,s( x ) dW(x). 

JR d 

By definition £f(r) = (/,1e t )- Therefore, 



(27) 



( T ;6)-l f (T)\ 2 <2\(f,Vr,5-lE T }\ 2 + 2 

<ci{d)M 2 5 2 + 2\X(u,r)\ 2 , 



e / i> T Ax)dW{x) 
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where we have taken into account the definition of <fi T .s(') and (12); X(u,r) 
is defined in (22). In view of Lemma 2, we have, for any t > 2ax, 



E sup \X(u,r)\ 2 = [ 2tp( sup \X(u, r)\ > t X dt 
re[o,i] Jo Ue[o,i] J 

/•oo 

< 2? + c 3 y_ iV^ 2 exp{-t 2 /(2o-|)} (it 



= 2? + C4<7x / v^sexpl— s} ds 

JP/2a 2 x 

< 2? + c 4 a x exp{-P/(4a x )}, 

where the last step uses Jensen's inequality. Minimizing with respect to i, 
we obtain Esup rg [ 0) i] \X(u,r)\ 2 < c^a x ln(l/cix), and the statement of the 
theorem follows from (23). □ 

Proof of Theorem 2. Throughout the proof we fix a direction u £ 
S^" 1 and therefore, for brevity, omit the argument u. With slight abuse of 
notation, we write h for h(u), if(r) for £f(u,r) and so on. We denote by 
c\ , C2 , . . . positive constants depending on d, and A only. 
For the estimate h*, we have 



E\h* - h\ 2 = E[\h* - h\ 2 l{h - A < h* < h}] 
(28) + E[\h, -h\ 2 l{h, <h- A}] 

1 ' +E[\K-h\ 2 l{K>h}] 

=:h+h + h, 

where A is defined in (17). We bound I\, 1% and I3 separately. 
In view of (17), we have 

\K - h\ 2 \{h -A<h*<h}< Q~ 2 / a \£ f (h*)\ 2/a 

< Q~ 2/a {\if(k) - l(h*)\ + \l(K)\} 2/a 

< ciQ- 2 / a {\£ f (K) - l(h*)\ 2/a + \l(k)\ 2/a } 

< Cl Q 2 / Q ( sup \l( r )-£ f (r)\ 2 / a + e 2 M, 
Lre[/i-A,ft] J 

(29) 

where we have used the fact that £*(/»*) = 9* with probability 1. Taking the 
expectation and using Theorem 1, we obtain 

< / I r\ 2/(a(/3+l/2)) n 

/i<c 2 Q" 2 / Q Ul^/ 2 y In-) +9 2 J a l 

For the second term I2, we have 

h = E[|/i* -h\ 2 l{K <h- A}] 
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< F{h* <h- A} 

< F{h(h-A)<6*} 

= F{£f(h-A)+i*(h-A)-£ f (h-A)<0*} 

< F{\£ f (h - A) - e*(h - A)| > £ f (h - A) - 0*} 

< F{\X(u,h-A)\ >£ f (h-A)/2} 

< ) P{|W(0,l)|>QA a ( e |^ Tj i.||)- 1 /2} 

( <c 3 exp{-C4 e - 2 /^ +1 / 2 )}, 

where (a) follows from definition of X(u,r), the choice of <5* and the fact 
that e is sufficiently small, and (b) follows from Lemma 1. Thus, I2 = o{I\) 
as e — > 0. 

It remains to bound the third error term 1%. We have 
I 3 = E[|/t* - h\ 2 l{K > h}} < F{K > h} 

= p{ sup 4(r)>0*}. 
Ue[ft,i] J 

Observe that 

4 (r) = (A A ,Kf)+e f A A (x) dW(x) 

= X(u,r) for r G (/i, 1], 

because the support of <p T> s t is E T , and E T n G = for r £ (/i, 1]. Note 
that, for e small enough, 0* > c^ox y/ln(l/ax)', see (23). Therefore, choosing 
9 = 9*, we have, by Lemma 2, 

/ 3 <p( sup \X(u,r)\ 

I rG [ft, 1] J 

< c^lnil/ax^axr- 1 < c 8 9 2 J a , 

where the last inequality is obtained by the choice of C3 in (18). Combining 
bounds for Ji, I2 and ^3, we complete the proof. □ 

Proof of Theorem 3. First we note that 

A 00 (G'*,G)< sup \h*(u) — h(u)\. 
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The remainder of the proof is identical to the proof of Theorem 2 with the 
following minor modifications. When bounding I\ in (28), we take the addi- 
tional supremum over u G S^ -1 . Therefore, the expectation Esup TgT |^*(t) — 
^/(t)| 2 should be bounded. This is done exactly as in the proof of Theorem 1, 
with the only difference being that the stochastic error of the estimator is 
now considered as the absolute value of the Gaussian process with the in- 
dex set T; then part (ii) of Lemma 1 is used. We note that in this case the 
constant C3 in (18) should be taken to be larger than in Theorem 1. In the 
same manner, when bounding J3, the statement (ii) of Lemma 1 is used. □ 

Acknowledgment. The authors would like to thank the referees for their 
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